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1.  Introduction 


This  report  presents  the  results  of  a  study  of  a  large  suite  of  T-phase  signals  from  submarine  earthquakes 
in  the  Pacific,  recorded  on  hydrophones  located  near  Pt.  Sur  and  Wake  Island.  The  records  were  obtained 
via  online  requests  from  the  Center  for  Monitoring  Research  (CMR)  for  an  8-10  minute  band  about 
each  expected  T-phase  arrival.  Event  parameters  were  obtained  from  the  near-real-time  Earthquake 
Bulletin,  provided  by  the  National  Earthquake  Information  Service  (NEIC).  The  events  cover  a  wide 
range  of  locations,  providing  good  azimuthal  coverage  from  source  to  receiver,  as  well  as  a  wide  range 
of  source  depths. 

The  motivations  for  this  study  were: 

1 .  to  determine  physical  mechanisms  by  which  energy  from  sources  at  various  depths  couples  into  the 
ocean  sound  channel  (or  SOFAR  channel),  and 

2.  to  determine  the  relative  importance  of  source  and  path  effects  in  detecting  T-phases. 

Since  volume  attenuation  of  acoustic  energy  is  negligible  at  frequencies  of  10-100  Hz  (Urick,  1979), 
the  sound  channel  is  an  extremely  efficient  waveguide  for  acoustic  propagation.  The  excellent  acoustic 
propagation  characteristics  of  the  sound  channel  allow  one  to  record  small  underwater  earthquakes  and 
nuclear  blasts  on  hydrophones  at  distances  of  thousands  of  kilometers  from  the  source.  A  number  of 
studies  (eg.  Duennebier,  1968;  Northrup,  1968, 1974;  Shurbet  and  Ewing,  1957;  Walker  et.aL,  1992; 
Walker  and  Bernard,  1993)  have  shown  that  T-phases  from  suboceanic  earthquakes  vary  considerably 
from  region  to  region,  and  possibly  with  magnitude  or  earthquake  source  mechanism.  Since  monitoring 
a  CTBT  requires  the  ability  to  distinguish  between  oceanic  nuclear  detonations  and  earthquakes,  it  is 
important  to  be  able  to  distinguish  between  T-phases  resulting  from  these  sources.  Toward  this  end,  we 
seek  to  understand  the  coupling  of  seismic  energy  to  acoustic  energy  in  the  sound  channel  and  the 
propagation  of  acoustic  energy  from  source  region  to  receiver. 

Our  smdy  is  composed  of  3  parts.  In  the  next  section,  we  present  a  study  of  T-phase  amplitudes  vs. 
earthquake  magnitude  and  event  location.  Since  instrument  calibration  information  was  not  available, 
an  amplitude  measure  was  devised  based  on  the  signal- to-noise  levels  at  each  station.  A  comparison  of 
T-phase  source  locations  determined  by  T-phase  arrivals  at  the  Wake  and  Pt.  Sur  hydrophones  with  the 
NEIC-predicted  epicentral  locations  is  presented  in  Section  3,  along  with  a  discussion  of  sytematic 
differences  between  these  sets  of  locations.  In  Section  4,  a  method  of  modeling  T-phase  coda  for  sub- 
oceanic  earthquakes  is  presented  which  is  based  on  the  assumption  that  the  coda  may  be  estimated  as 
the  summation  of  the  seafloor  displacements  in  the  source  region. 

2.  Dependence  of  T-phase  amplitude  on  source  strength  and  event  location 

2.1  A  T-phase  amplitude  measure 

Since  acoustic  energy  loss  in  the  ocean’s  sound  channel  is  negligible,  the  energy  recorded  at  hydrophones 
at  distances  of  several  thousand  kilometers  from  the  source  is  expected  to  be  directly  proportional  to 
the  seismic  energy  coupled  to  the  sound  channel  for  paths  which  do  not  intersect  ridges  or  island 
chains. 
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Sonograms  of  T-phases  generated  by  suboceanic  earthquakes  indicate  that  most  of  the  acoustic  energy 
IS  concentrated  between  2  and  20Hz.  It  seems  probable  that  the  high  magnitudes  in  this  frequency 
range  are  partly  due  to  the  instrument  response  since,  for  instance,  one  might  expect  to  see  microseismic 
energy  at  very  low  frequencies.  Since  we  do  not  have  instrument  calibration  information,  we  cannot 
convert  these  data  into  absolute  pressure  magnitudes.  Instead  we  introduce  an  amplitude  measune  of 
the  T-phase  which  depends  on  observed  SNR  ratios. 

The  method  by  which  the  T-phase  amplitude  is  defined  is  shown  in  Figure  1  for  several  frequency 
bands  for  a  magnitude  5.6  Aleutian  Island  event  (located  at  51.35°N,  183.49°  E).  Band-passed  filtered 
versions  of  the  data  as  recorded  at  the  Pt  Sur  hydrophone  are  shown  for  5  frequency  ranges.  For  each 
filtered  trace,  the  signal  variance  was  determined  over  5  second  band,  with  steps  of  1  second  between 
bands.  A  ratio  R,  given  by  the  maximum  variance  divided  by  the  minimum  variance,  was  computed  for 
each  filtered  trace.  The  signal-to-noise  ratio  SNR  is  given  by  R-1.  Scaled  SNRs  are  shown  to  the  right 
for  each  frequency  band.  The  T-phase  amplitude,  which  is  defined  as  the  maximum  SNR  for  each 
frequency  band,  is  shown  for  each  trace.  For  this  event,  the  maximum  signal  to  noise  ratio  is  in  the  4- 


Figure.  1  The  T-phase  amplitude  is  defined  using  the  signal-to-noise  ratio  and  is 
sho\vn  here  for  several  frequency  bands.  The  filtered  trace  are  shown  to  the  left, 
the  signal  variance  as  a  function  of  time  is  shown  to  the  left.  T-phase  amplitudes 
are  given  for  each  band. 
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8  Hz  frequency  band.  For  most  events,  the  T-phase  signal  is  primarily  in  the  2  to  8  Hz  frequency  band, 
therefore  this  band  was  used  for  computation  of  T-phase  amplitudes  for  all  events  in  the  next  section. 

2.2  T-phase  amplitudes  observed  at  the  Pt.Sur  and  Wake  Island  hydrophones 
In  Figure  2,  events  corresponding  to  T-phases  observed  at  the  Pt  Sur  hydrophones  are  indicated  by 
circles  plotted  at  the  event  locations.  The  circle  size  is  an  increasing  function  of  the  observed  T-phase 
amplitude.  Hydrophone  locations  are  shown  as  asterisks,  events  with  no  observable  T-phase  energy  are 
shown  as  x’es.  The  corresponding  diagram  for  T-phases  observed  at  Wake  is  shown  in  Figure  3.  Note 


T  phase  amplitudes  at  Pt.  Sur 


Figure.  2  T-phase  amplitudes  observed  at  Pt.  Sur,  shown  as  a  function  of 
earthquake  location.  The  circle  size  is  proportional  to  the  square  root  of  the 
amplitude.  X’es  indicate  earthquakes  for  which  no  T-phase  was  discernible.  Asterisks 
indicate  hydrophone  location. 


that  a  large  number  of  events,  having  varying  source  magnitudes,  depths,  source  mechanisms,  and 
transmission  path  effects  have  been  included  in  the  above  analysis. 

In  Figure  2,  the  geodesic  path  which  intersects  the  Aleutian  chain  is  shown  as  an  arc.  Calculations  of 
acoustic  transmission  along  geodesic  paths  indicate  that  events  to  the  northeast  of  this  line  lie  in  a 
shadow  region  with  respect  to  the  Pt  Sur  hydrophone,  i.e.,  acoustic  transmission  for  these  events  is 
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T  phase  amplitudes  at  Wake 


longitude 


Figure.  3  T-phase  amplitudes  observed  at  Wake  Island. 


blocked  by  the  Aleutian  chain.  As  shown,  there  is  little  or  no  discernible  T-phase  energy  at  the  Pt  Sur 
hydrophone  for  events  near  Kamchatka  and  in  the  Kuril  Island  chain.  These  events  do,  in  fact,  generate 
a  considerable  amount  of  T-phase  energy,  as  indicated  in  the  corresponding  figure  for  transmissions 
observed  on  the  Wake  hydrophone.  There  are  no  significant  bathymetric  features  blocking  acoustic 
uansmission  to  Wake  Island  for  these  events,  and  T-phase  energy  is  readily  discernible  at  the  Wake 
hydrophone.  Thus,  T-phase  amplitudes  and  coda  for  this  region  are  strongly  dependent  upon  path 
effects,  rather  than  source  effects.  For  events  which  are  readily  observable  at  both  hydrophones  lower 
noise  levels  at  Wake  generally  result  in  greater  T-phase  amplitude  measures  there  than  at  Pt.  Sur. 

2.3  Dependence  of  T-phase  amplitudes  on  source  magnitude 

In  order  to  determine  how  the  T-phase  amplitude  scales  with  event  magnitude,  we  examine  a  cluster  of 
events  which  occured  near  the  Aleutian  Islands  in  the  latitude  range  51“  to  52“  North,  and  longitude 
range  from  181°  to  184“  East.  By  examining  a  cluster  of  events  we  eliminate  effects  due  to  varying 
distances,  earthquake  types,  event  depths  and  varying  amounts  of  transmission  blockage.  The  event 
locations  were  taken  from  the  NEIC  and  all  locations  accuracies  were  rated  A  or  B  All  events  are 
likely  associated  with  the  subduction  of  the  Pacific  Plate  at  the  Aleutian  trench  and  occur  at  shallow 
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depths  (between  28  and  60  km).  In  all,  this  dataset  is  composed  of  33  events.  There  is  little  transmission 
blockage  along  the  geodesic  paths  from  this  region  to  PL  Sur  or  to  Wake. 


The  Log,,  (T-phase  amplitude)  values  as  observed  at  the  Pt  Sur  hydrophone  are  plotted  as  a  function 
of  magnitude  in  Figure  4  for  all  five  frequency  bands.  The  amplitudes  are  an  increasing  function  of 
magnitude  as  expected,  but  there  remains  a  considerable  amount  of  scatter.  The  2  to  4  Hz  frequency 
band  shows  the  greatest  increase  in  T-phase  amplitude  with  event  magnitude,  and  the  least  scatter. 
However,  the  T-phase  amplitudes  in  the  4  to  8  Hz  band  are  slightly  larger  on  average.  At  the  very 


T  phase  amplitude  vs.  magnitude  for  Aleutian  Island  cluster 


Figure.  4  T-phase  amplitudes  observed  at  Pt.  Sur,  as  a  function  of  magnitude  for 
a  cluster  of  earthquakes  near  the  Aleutian  Islands.  The  greatest  correlation  between 
T-phase  amplitude  and  magnitude  is  in  the  2-4  Hz,  and  4-8  Hz  frequency  bands. 


highest  frequencies  examined,  16  to  32  Hz,  there  is  little  correlation  between  magnitude  and  T-phase 
amplitude. 

The  corresponding  plot  for  observations  at  the  Wake  hydrophone  is  shown  in  Figure  5.  Also  plotted  are 
the  best  fit  lines  for  each  frequency  range.  Note  that  there  is  a  large  degree  of  scatter  about  the  best  line 
fit  and,  as  for  Pt  Sur,  the  T-phase  amplitudes  are  greatest  in  the  2  to  8  Hz  frequency  range.  However,  the 
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T  phase  amplitude  vs.  magnitude  at  Aleutian  Island  cluster,  observed  at  Wake 


Figure.  5  As  in  Figure.  4,  but  for  T-waves  observed  at  Wake.  Best  fit  lines  for  each 
frequency  band  are  shown. 


slope  of  each  line  decreases  with  increasing  frequency.  In  general,  T-phase  amplitudes  are  somewhat 
higher  at  Wake  than  at  Pt  Sur,  probably  due  to  lower  noise  levels  at  Wake. 

T-phases  are  shown  in  Figure  6  as  a  function  of  event  magnitude  for  a  sub-cluster  of  the  events,  (from 
183°  to  184°  East).  Labels  to  the  right  of  the  traces  show  the  distance  in  degrees  from  the  event  to  the 
Pt  Sur  hydrophone.  As  indicated,  the  T-phase  signatures  are  fairly  uniform  as  a  function  of  magnitude 
within  a  given  geographic  area, although  the  duration  and  amplitude  of  the  T-phase  appears  to  scale 
with  magnitude  within  this  limited  region.  The  scatter  in  amplitude  and  duration  is  significant,  however. 

3.  T-phase  source  locations 

One  purpose  of  the  hydroacoustic  monitoring  system  for  the  CTBTis  to  detect  and  locate  both  underwater 
events  and  suboceanic  events.  Although  observations  from  at  least  three  stations  are  usually  required 
to  accurately  locate  an  event,  only  two  hydrocoustic  stations,  at  Wake  Island  and  Pt  Sur,  are  in  place 
for  detection  of  hydroacoustic  events  in  the  Pacific.  Thus,  event  source  location  in  the  Pacific  will 
generally  rely  on  the  analysis  of  a  combination  of  seismic  and  hydroacoustic  data  types.  However,  it 
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T  phase  signature  at  various  magnitudes 


321  0123456 


reduced  travel  time  (min) ;  Tvel  =  1 .5km/sec 

Figure.  6  T-phases  for  an  Aleutian  Island  cluster  filtered  to  2-8  Hz.  Note  that, 
although  amplitude  and  duration  scale  with  magnitude,  the  character  of  the  T- 
phase  coda  is  similar  for  earthquakes  within  this  cluster. 

has  been  noted  by  numerous  authors  {e.g.  Johnson,  1966)  that  the  seismic/T-phase  coupling  may  occur 
at  locations  other  than  at  the  earthquake  epicenters.  Therefore,  the  success  of  an  integrated  seismic/ 
hydroacoustic  analysis  will  depend  on  whether  the  locus  of  T-phase  excitation  is  coincident  with  the 
locus  of  P  and  S-wave  excitation,  i.e.  at  the  earthquake  epicenter.  In  this  section,  we  determine  T-phase 
source  locations  using  hydroacoustic  data  from  Pt  Sur  and  Wake,  and  compare  them  with  source  locations 
determined  by  the  NEIC.  This  wiU  indicate  the  magnitude  of  location  error  which  results  from  assuming 
that  the  T-phase  source  is  coincident  with  the  earthquake  epicentral  location. 

Johnson  (1966)  published  a  method  for  “routinely”  determining  T-phase  source  locations,  but  noted 
that  these  locations  were  not  necessarily  coincident  with  earthquake  epicentral  locations.  In  this  method, 
the  T-phase  onset  time  was  taken  as  the  largest  T-phase  amplitude  in  the  frequency  range  from  1-  15Hz. 
Application  of  Johnson’s  method  to  a  large  series  of  earthquakes  with  known  epicentral  locations 
occurring  over  a  fairly  uniform  area  showed  that  the  T-phase  energy  appeared  to  be  generated  at  several 
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distinct  regions  (Johnson  and  Norris,  1968),  suggesting  that  T-phases  were  being  generated  only  at 
isolated  locations  along  the  ocean  floor.  Furthermore,  it  has  been  shown  that  seamounts  and  islands  act 
as  secondary  generators  of  low  frequency  T-phase  energy  (Walker  et  al.,  1992;  Northrup,  1973,1974). 
Thus,  accurate  source  location  becomes  complicated  in  locations  where  phases  associated  with  secondary 
sources  are  of  the  same  order  of  magnitude  as  the  primary  T-phase  radiator. 

A  suite  of  records  from  54  events  for  which  hydroacoustic  data  were  available  from  stations  at  both  Pt. 
Sur  and  Wake  Island  was  assembled.  We  computed  envelopes  of  the  T-phase  coda  for  several  frequency 
ranges.  We  noted  that  for  several  events,  T-phase  onset  times,  as  defined  by  the  time  of  the  largest  T- 
phase  amplitude,  were  strongly  dependent  on  the  frequency  range  examined.  This  occurred  where 
secondary  sources  generated  larger  low  frequency  T-phases  than  the  primary,  i.e.  near-source,  radiator. 
To  resolve  the  ambiguity  in  T-phase  onset,  we  chose  the  maximum  amplitude  arrival  in  the  10-20Hz 
frequency  range.  Onset  time  picks  are  accurate  to  2  seconds. 

To  determine  travel  times  for  arbitrary  source  locations  to  the  receivers  at  Wake  and  Pt  Sur,  we  integrated 
SOFAR  channel  velocities  along  geodesic  paths.  Sound  speeds  were  computed  from  seasonal  salinity 
and  temperature  databases,  provided  by  Farrell  and  LePage  (1996).  We  used  the  annual  averages  of  the 
SOFAR  channel  velocities;  negligible  differences  were  noted  between  travel  times  computed  at  different 
seasons.  Note  that  3-D  acoustic  effects,  i.e.  the  horizontal  refraction  of  acoustic  energy  near  islands 
and  seamounts,  were  not  taken  into  account 


Travel  times  to  Wake  and  Pt  Sur  are  shown  in  the  top  two  panels  of  Figure  7.  Shadow  zones  regions  are 
roughly  indicated  by  the  white  patches;  these  are  regions  where  the  bathymetry  along  the  geodesic 
path  intersects  the  ocean  surface.  Travel  time  differences  are  shown  in  the  bottom  panel  of  Figure  7. 
Note  that,  without  supplemental  seismic  or  T-phase  information  (from  island  T-phase  stations,  for 
instance),  the  difference  in  T-phase  onset  times  at  Pt  Sur  and  Wake  are  the  only  information  available 
to  determine  source  location.  This  will  be  the  case  for  most  small  oceanic  and  sub-oceanic  sources 
because  T-phases  can  generally  be  more  reliably  detected  than  seismic  body  waves  (Brocher,  1983; 
Shurbet,  1962),  due  to  the  extremely  efficient  propagation  of  acoustic  energy  through  the  ocean.  In 
this  case,  note  that  the  longitude  of  the  T-phase  source  is  much  more  clearly  resolved  than  the  latitude 
over  most  of  the  Pacific.  Resolution  is  even  worse  for  regions  to  the  southeast  of  Wake  Island. 

T-phase  source  locations  were  determined  using  the  travel  times  as  given  above,  as  well  as  the  given 
NEIC  origin  times  and  T-phase  onset  times  defined  above.  Travel  times  through  the  crust  were  neglected 
in  determining  source  locations.  NEIC  epicenters  and  T-phase  source  locations  are  compared  in  Table 
1,  and  plotted  in  Figure  8.  For  several  events  near  coasts,  no  accurate  T-phase  source  location  could  be 
determined.  Two  possible  causes  for  this  are  that  either  (1)  horizontal  refraction  must  be  accounted  for 
near  the  coasts,  or  that  (2)  peaks  corresponding  to  a  near-source  radiator  was  picked  for  one  station, 
whereas  a  peak  corresponding  to  a  secondary  source  was  picked  for  the  other  station.  Several  NEIC 
events,  indicated  by  asterisks,  were  not  located  using  T-phases  as  they  were  not  observable  at  one  of 
the  stations.  Figure  9  shows  a  blowup  of  the  epicentral  locations  and  T-phase  source  locations  in  the 
Aleutian  Islands,  superimposed  on  a  bath)metric  map.  Note  that  T-phase  source  locations  are  consistently 
biased  to  the  north  of  the  NEIC  locations,  in  regions  of  shallower  water. 
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Figure.  7  Travel  times  to  Pt.  Sur 
and  Wake  are  shown  in  the  top 
two  panels.  The  travel  time 
difference  is  shown  in  the  bottom 
panel. 
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Figure.  8  NEIC  locations  (x’es)  compared  to  apparent  T-phase  source  location 
(circles)  for  the  North  Pacific. 


map  or  bathymetry.  The  apparent  T-phase  sources  are  biased  toward 
shallower  regions  than  the  epicentral  locations. 


Table  1:  NEIC  and  T-phase  locations 


Julian 

origin  time 

epicentral 

epicentral  depth  mag 

quality  T-source  T-source 

day  1996  (hr:min:sec) 

latitude(°N)  longC 

=E)  (km) 

latitude  long(°E) 

156 

0:37:18 

46.84 

153.70 

33.0 

4.6 

A 

no  signal  at  Pt.  Sur 

156 

3:30:47 

52.90 

192.30 

33.0 

4.2 

A 

53.2 

192.8 

163 

0:18:46 

51.22 

183.37 

33.0 

4.6 

B 

51.95 

182.85 

163 

2:14:57 

51.02 

183.36 

33.0 

4.5 

B 

52.05 

182.95 

163 

3:44:40 

51.31 

184.08 

33.0 

4.9 

B 

52 

183.2 

163 

5:44:02 

52.14 

181.65 

33.0 

4.5 

B 

51.65 

181.85 

163 

10:40:08 

51.30 

183.76 

33.0 

5.7 

A 

52.25 

183.2 

163 

13:03:03 

51.35 

183.49 

33.0 

5.6 

A 

51.85 

183.55 

163 

18:22:55 

12.71 

125.00 

33.0 

7.0 

B 

10.80 

126.3 

165 

2:47:29 

51.17 

182.44 

33.0 

4.4 

B 

51.95 

183.2 

168 

15:01:29 

51.31 

183.10 

33.0 

5.1 

A 

51.85 

183.2 

168 

19:31:16 

51.42 

181.45 

30.2 

5.3 

A 

51.6 

181.65 

169 

4:15:54 

5.11 

277.63 

10.0 

4.8 

B 

no  signal  at  Wake 

169 

16:00:47 

51.33 

181.73 

33.0 

5.0 

A 

51.7 

181.8 

169 

21:55:35 

13.81 

143.92 

33.0 

5.4 

B 

9.0 

147.4 

170 

16:06:24 

51.26 

182.72 

33.0 

4.8 

B 

52.05 

182.85 

171 

3:20:26 

4.90 

277.55 

33.0 

5.0 

A 

no  signal  at  Pt.  Sur 

171 

8:33:15 

51.19 

182.41 

33.0 

4.5 

B 

no  signal  at  Pt.  Sur 

171 

13:03:02 

43.49 

147.57 

33.0 

4.9 

B 

43.6 

147.6 

173 

13:57:11 

51.77 

159.01 

33.0 

6.6 

B 

51.4 

159.9 

174 

3:43:58 

51.60 

159.46 

33.0 

4.9 

B 

no  signal  at  Pt.  Sur 

178 

3:22:03 

27.74 

139.83 

468.0 

5.4 

A 

no  signal  at  Wake 

178 

12:35:55 

52.70 

192.62 

33.0 

4.4 

B 

53.25 

192.7 

182 

0:52:06 

51.34 

182.67 

33.0 

5.1 

B 

51.7 

183.2 

182 

9:51:30 

22.20 

143.05 

268.9 

4.8 

B 

23.5 

142.5 

182 

11:32:35 

51.02 

159.71 

33.0 

5.9 

B 

51.4 

161 

182 

23:42:15 

53.81 

169.71 

33.0 

5.2 

A 

no  signal  at  Pt. 

Sur 

187 

5:15:05 

51.32 

183.11 

33.0 

4.5 

B 

51.8 

183.35 

188 

5:03:33 

3.94 

128.31 

33.0 

5.5 

B 

no  signal  at  Pt. 

Sur 

188 

11:56:44 

15.72 

147.52 

33.0 

5.7 

A 

13.8 

149.4 

188 

21:36:28 

22.06 

142.80 

240.5 

5.7 

A 

13.8 

149.4 

189 

4:19:33 

53.14 

193.87 

53.2 

4.7 

B 

53.75 

193.85 

191 

7:56:43 

52.15 

189.94 

33.0 

4.2 

B 

53.45 

188.75 

191 

12:02:47 

51.90 

151.55 

500.0 

5.1 

A 

48.8 

150.4 

192 

0:04:21 

51.15 

183.40 

33.0 

4.7 

B 

52.05 

182.8 

192 

0:51:05 

50.91 

187.23 

33.0 

4.4 

B 

no  signal  at  Pt. 

Sur 

192 

5:48:19 

51.99 

188.86 

33.0 

5.4 

A 

52.65 

189.20 

192 

8:33:09 

51.76 

188.97 

33.0 

4.2 

B 

52.60 

189.30 

195 

15:10:30 

51.30 

182.10 

33.0 

5.3 

A 

51.70 

182 

198 

2:32:47 

45.90 

149.86 

33.0 

5.0 

B 

no  signal  at  Pt. 

Sur 

198 

3:48:27 

56.00 

165.00 

33.0 

6.3 

A 

no  signal  at  Pt. 

Sur 

198 

6:12:13 

56.06 

164.57 

33.0 

4.9 

B 

no  signal  at  Pt. 

Sur 

198 

9:00:41 

55.94 

164.55 

33.0 

4.8 

B 

no  signal  at  Pt. 

Sur 

200 

8:55:21 

51.13 

183.67 

33.0 

4.5 

B 

no  signal  at  Pt. 

Sur 

200 

9:26:53 

51.95 

179.55 

150.0 

4.6 

B 

no  signal  at  Pt. 

Sur 

204 

14:19:35 

0.85 

120.15 

33.0 

6.9 

B 

1.3 

120.9 

205 

13:24:56 

18.58 

204.35 

10.0 

4.7 

B 

17.9 

205.5 

206 

0:49:37 

33.30 

137.15 

33.0 

5.1 

B 

no  signal  at  Pt. 

Sur 

206 

20:15:44 

41.77 

234.01 

10.0 

5.5 

A 

42.3 

205.5 

206 

23:53:22 

51.64 

182.80 

59.4 

5.2 

A 

51.70 

183.15 

207 

9:30:20 

50.77 

177.40 

33.0 

4.6 

B 

52.4 

176.8 

210 

9:30:21 

18.83 

204.68 

10.0 

5.0 

B 

18.00 

205.5 

211 

20:20:55 

25.29 

122.26 

73.3 

5.3 

B 

23.7 

123.8 
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4.  Modeling  T-phase  coda 


In  this  section,  we  demonstrate  an  approximate  method  of  simulating  the  T-phase  coda.  This  may  lead 
to  methods  of  T-phase  source  location  using  only  one  or  two  receivers.  It  has  been  observed  by  Walker 
e£.a/.(1992)  that  the  character  of  the  T-phase  coda  is  strongly  dependent  upon  the  bathymetry  of  the 
ocean  floor  in  the  region  in  which  seismic  to  acoustic  coupling  takes  place.  Thus,  earthquakes  from  the 
same  geographic  area  generate  similar  T-phase  coda.  This  is  demonstrated  in  Figures  10  and  11,  which 
show  sonograms  from  a  cluster  of  events  in  the  Aleutian  Islands.  The  epicentral  locations  are  listed  in 
Table  2,  in  order  from  west-to-east.  Note  that  T-phases  observed  at  Pt.  Sur  (Figure  10)  have  a  relatively 
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Figure.  10  Sonograms  for  T-waves  observed  for  Pt.  Sur,  for  earthquakes  in  a 
cluster  from  51°to  52°N,  ISPto  185°E.  Epicentral  locations  are  listed  in  Table  2. 
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slow  rise  time,  and  are  quite  symmetric  about  the  defined  onset  time.  The  T-phases  at  Wake  (Figure  1 1 ) 
are  much  less  symmetric,  have  a  sharper  rise  time  and  a  longer  duration  than  at  Pt  Sur. 

Table  2:  NEIC  epicentral  locations  for  Aleutian  Island  cluster 


latitude 

longitude 

labels  in  figures 

(°N) 

(°E) 

10  and  11 

51.42 

181.4 

a 

51.33 

181.7 

b 

51.30 

182.1 

c 

51.34 

182.7 

d 

51.64 

182.8 

e 

51.31 

183.1 

f 

51.02 

183.4 

g 

51.35 

183.5 

h 

51.30 

183.8 

i 

51.31 

184.1 

j 

This  observation,  that  the  T-phase  coda  depends  upon  source  bathymetry  rather  than  path  effects  is 
further  supported  by  Figure  12.  T-phases  observed  at  Wake  for  events  at  several  azimuths  are  shown 
here  as  a  function  of  reduced  travel  time.  The  reducing  velocity  was  1.5km/sec.  The  traces  have  been 
band-passed  to  2-8Hz.  All  events  are  shallow  and  are  within  30  to  40°  from  Wake.  Event  magnitudes 
are  labelled  to  the  right  of  each  trace.  Note  that  while  events  within  a  cluster  generate  T-phases  with 
similar  characteristics,  these  characteristics  vary  from  cluster  to  cluster.  Both  the  duration  of  the  T- 
phase  and  its  onset  time  (with  respect  to  the  expected  arrival  at  t=0)  vary  considerably  with  azimuth 
For  instance,  the  onset  time  for  the  event  to  the  east  of  Wake  (at  27.5°  N,  128.3°  E)  is  considerably 
advanced  with  respect  to  the  expected  onset,  suggesting  that  conversion  from  seismic  to  acoustic  energy 
occurs  at  a  location  nearer  to  the  hydrophone.  These  variations  in  T-phase  characteristics  with  azimuth 
may  be  useful  for  determining  earthquake  locations  using  a  small  number  of  receivers. 

T-phase  envelopes  computed  over  2  Hz  frequency  bands  from  the  Wake  and  Pt  Sur  records  for  an  event 
at  51.99N,  188.86E  are  shown  in  the  top  and  bottom  panels  of  Figure  13,  respectively.  This  event  is 
approximately  equidistant  from  the  Wake  and  Pt.  Sur  hydrophones.  Note  that  there  are  two  main  peaks 
in  the  T-phase  amplitudes  at  each  station.  This  does  not  appear  to  be  a  source  effect  as  the  larger  peak 
follows  the  smaller  peak  at  Wake  and  precedes  it  at  Pt  Sur.  We  hypothesize  that  these  peaks  are  due  to 
seismic/acoustic  coupling  occuring  over  a  wide  area  of  the  seafloor.  Bathymetry  for  this  region  is 
shown  in  Figure  14,  with  superimposed  contours  of  the  travel  times  to  Pt  Sur  (solid  line)  and  to  Wake 
(dashed  line).  Travel  times  from  the  hypocenter  are  given  by  the  summation  of  T-phase  travel  times  and 
the  crustal  travel  time  (assuming  a  crustal  velocity  of  8km/sec).  The  epicenter  is  shown  by  the  black 
circle.  Note  that  the  larger  T-phase  peaks  on  the  Pt  Sur  and  Wake  records  corresponds  to  a  T-phase 
source  location  at  approximately  52.5N,  189.5E,  i.e.  at  a  bathymetric  high  to  the  northeast  of  the 
epicenter;  the  smaller  peak  corresponds  to  a  T-phase  source  at  about  52.2N,  188E,  i.e.  at  a  bathymetric 
high  to  the  northwest  of  the  epicenter. 

If  the  T-phase  source  location  were  coincident  with  the  epicentral  location,  the  expected  travel  time  to 
Wake  would  be  2790  sec,  and  to  Pt  Sur  would  be  2840  seconds.  As  indicated  in  Figure  13,  the  T-phase 
onset  at  Wake  is  delayed  with  respect  to  the  expected  travel  time.  Thus,  the  onset  time  is  consistent 
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Figure.  1 2  T-phase  coda  for  a  variety  of  earthquakes  at  epicentral  distances  of  30°  to 
40°  from  Wake.  Note  that,  although  the  character  of  T-phase  coda  is  similar  within 
regions,  it  varies  considerably  from  one  region  to  another,  as  indicated. 
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Figure.  14  Bathymetry  map,  with  superimposed  travel  times  to  Pt.  Sur  (solid 
lines)  and  to  Wake  (dashed  lines).  Travel  times  are  computed  for  the  event  at 
51.99°N,  188.86°E,  with  a  depth  of  33km,  assuming  a  crustal  velocity  of  8km/sec 
and  SOFAR  channel  velocities  along  geodesic  paths. 


with  seismic  energy  travelling  through  the  crust  to  a  shallower  region  slightly  north  of  the  epicenter, 
then  being  converted  to  acoustic  energy. 

Several  authors  have  theorized  (Tolstoy  and  Ewing,  1950;  Johnson  et  al.,  1963)  that  transformation  of 
P  to  T  waves  takes  place  along  the  continental  slope.  Several  studies  (Okal  and  Talandier,  1986; 
Johnson  et.al,  1963;  Ewing  et.al,  1950)  have  argued  that  the  trapping  of  seismic  energy  into  the  SOFAR 
channel  requires  a  strongly  dipping  ocean  floor.  We  hypothesize  that  seismic/acoustic  coupling  occurs 
over  a  distributed  region  of  the  ocean  floor,  with  conversion  to  low  order  modes,  which  propagate  more 
efficiently  over  large  distances,  occurring  in  shallow  water.  Our  physical  model  is  that  of  seismic  to 
acoustic  coupling  occurring  mainly  by  scattering  at  the  ocean  floor,  with  downward  propagation  ii.e. 
mode  coupling  from  higher  to  lower  modes)  playing  a  negligible  role. 

We  choose  the  following  form  for  the  amplitude  of  T-phase  energy  originating  at  each  point  in  the 
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region  of  the  earthquake  epicenter: 


p(r,z)  =  g(z,f)  T(r)/d(r) 

where  p(r,z)  is  the  pressure  amplitude  at  the  seafloor,  d(r)  is  the  crustal  distance  from  the  epicenter  to 
the  ocean  bottom,  T(r)  is  the  transmission  factor  from  source  to  receiver  (where  T(r)  =  0  indicates  that 
no  energy  from  location 

r  propagates  to  the  receiver),  and  the  depth  factor  g(z,f)  is  a  function  of  ocean  bottom  depth  z  and 
frequency  f.  In  other  words,  we  compute  the  motion  at  each  point  along  the  seafloor  and  then  consider 
these  points  as  secondary  sources. 

Various  choices  may  be  made  for  the  form  of  g(z,f).  Since  higher  order  modes  will  generally  be  attenuated 
through  bottom  interaction,  we  take  g(z,f)  to  be  the  pressure  amplitude  of  mode  1  at  the  ocean  bottom. 
A  more  complete  description  of  g(z,f)  T(r)  would  include  a  summation  over  modes,  with  a  separate 
value  of  depth  factor  g  and  transmission  factor  T  for  each  mode.  However,  since  higher  modes  will 
undergo  greater  bottom  interaction,  and  hence  greater  attenuation,  we  neglect  all  modes  above  mode  1 
for  this  computation.  Assuming  a  seafloor  with  a  velocity  profile  given  by 
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the  mode  1  pressure  amplitude  g(z,f)  at  the  ocean  bottom  is  as  shown  in  Figure  15.  This  profile  was 
computed  for  a  seawater  velocity  profile  at  SON,  180E.  The  pressure  dropoff  with  depth  is  more 
gradual  for  more  southerly  velocity  profiles. 

Transmission  factors  for  Wake  and  Pt  Sur  are  shown  in  the  top  panels  of  Figure  16  for  the  event  at 
51.99N,  188. 86E.  White  areas  indicate  shadow  regions,  black  areas  indicate  regions  in  which  energy 
may  propagate  to  the  receiver.  The  pressure  amplitudes  at  the  seafloor,  given  by  g(z,f)/  d(r)  are  shown 
in  Figure  17,  with  darker  areas  corresponding  to  regions  of  higher  pressure  amplitudes.  Note  that  there 
are  2  larger  and  1  smaller  peak  in  these  pressure  amplitudes,  which  appear  to  nearly  coincide  with  the 
expected  locations  of  the  two  T-sources. 

To  compute  the  synthetic  T-phase  envelopes  we  sum  amplitudes  within  5  second  bins.  The  comparison 
of  real  and  synthetic  amplitudes  for  a  frequency  of  5  Hz  is  shown  in  Figure  17  at  both  Wake  and  Pt  Sur. 
Note  that  the  synthetic  envelope  (solid  line)  at  Wake  is  advanced  by  18  seconds  with  respect  to  the  real 
envelope,  but  the  overall  shape  of  the  envelopes  correspond  well.  Agreement  is  even  closer  for  the  Pt 
Sur  real  and  synthetic  envelopes.  Also  note  that,  in  computing  these  envelopes,  the  duration  and  spatial 
extent  of  the  source  were  not  taken  into  account. 

Similar  computations  for  sources  near  a  coastline  (not  illustrated  here)  indicated  the  source  location  of 
peak  energy  at  Wake  does  not  necessarily  coincide  with  that  of  Pt  Sur.  The  source  location  of  the 
maximum  T-phase  amplitude  at  each  hydrophone  depends  on  both  the  transmission  factor  {i.e.  blockage 
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pressure  amplitude  at  ocean  bottom  as  a  function  of  bottom  depth 


Figure.  1 5  Mode  1  pressure  amplitudes  at  the  ocean  bottom  for  f=5  and  10  Hz. 
This  profile  was  computed  for  a  seawater  velocity  profile  at  50°N,  180°,  and  a 
crustal  velocity  profile  as  given  in  the  text. 


from  source  to  receiver)  and  the  direction  of  propagation. 

The  comparisons  of  real  and  synthetic  T-phase  coda  amplitudes  illustrate  the  need  for  more  accurate 
determination  of  travel  times,  which  would  involve  the  use  of  mode  1  phase  speeds  instead  of  SOFAR 
channel  velocities,  and  also  the  effects  of  lateral  refraction.  The  difference  between  mode  1  phase 
speeds  and  the  SOFAR  channel  sound  speeds  is  shown  in  the  top  left  figure.  Note  that  there  is  an 
appreciable  difference  between  the  two  only  at  high  latitudes.  The  effect  of  lateral  refraction,  not 
shown  here,  would  be  to  bend  rays  further  to  the  south.  The  difference  between  mode  2  and  mode  1 
velocities  is  shown  in  Figure  18.  There  is  little  difference  between  the  2  speeds,  indicating  that  the 
signal  duration  can’t  be  accounted  for  by  mode  dispersion. 
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16  Transmission  factors  for  Wake  and  Pt.  Sur  are  shown  in  the  top  panels. 
V^ite  areas  indicate  shadow  zones,  which  correspond  to  transmission  factors  of  0 

The  bottom  profile  indicates  pressure  amplitude  at  the  seafloor  for  scattering  into 
mode  1.  ^ 
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normalized  amplitude  normalized  amplitude 


comparison  of  real  and  synthetic  envelopes  at  Wake  receiver 


Figure.  17  Comparison  of  real  (solid  line)  and  synthetic  (dashed  line)  amplitudes 
at  Wake,  top,  and  Pt.  Sur,  bottom.  Synthetic  envelopes  are  computed  at  5Hz;  real 
envelopes  are  computed  for  a  4-6  Hz  frequency  band. 
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longitude  (E) 


Figure.  18  A  comparison  of  mode  1  phase  velocities  and  SOFAR  channel  velocities 
is  shown  in  the  top  panel.  The  comparison  of  mode  1  and  2  phase  speeds  in 
shown  in  the  bottom  panel. 


5.  Conclusions 


The  conclusions  are  listed  in  point  form  below: 

•  T-phase  energy  drops  off  rapidly  at  frequencies  above  16  Hz. 

•  T-wave  coda  from  the  same  geographic  region  tend  to  be  very  similar.  Variations  in  T-phase  coda 
are  more  strongly  dependent  upon  source  region  than  on  source  mechanism,  depth,  or  magnitude. 
This  could  serve  to  determine  earthquake  epicentral  location  within  several  degrees. 

•  Although,  at  some  frequencies,  there  is  a  direct  scaling  between  T-phase  amplitude  and  source 
magnitude,  the  scatter  on  amplitude  vs.  magnitude  plots  is  too  great  to  use  amplitude  as  a  diagnostic 
measure  of  source  magnitude. 

•  T-phase  sources  in  the  Pacific  cannot  be  accurately  located  using  only  onset  times  at  the  Pt  Sur  and 
Wake  stations. 

•  There  is  a  systematic  difference  between  T-phase  source  location  and  earthquake  epicentral  location 
The  acoustic  energy  is  excited  in  shallow  areas.  This  implies  that  only  an  approximate  epicentral 
location  can  be  determined  from  a  joint  analysis  of  seismic  and  hydroacoustic  data. 

•  T-phase  source  location  tends  to  be  less  accurate  near  ocean/continent  boundaries.  This  may  be  due 
either  to  generation  of  T-phases  at  several  locations  along  the  coast,  with  different  amounts  of 
transmission  blockage  between  sources  and  receivers,  or  to  differences  in  propagation  direction 
from  source  to  receiver.  That  is,  an  elongated  region  of  high  seafloor  pressure  amplitudes  which  is 
perpendicular  to  the  propagation  direction  to  a  particular  hydrophone  may  show  up  as  a  maximum 
on  that  receiver,  but  show  up  as  coda. 

•  The  T-phase  onset  time  is  consistent  with  seismic/T-phase  conversion  at  relatively  nearby  shallow 
structures.  The  lengthened  arrivals  are  consistent  with  a  super-position  of  T-phases  generated  over 
shallow  regions  near  the  earthquake  epicenter. 

•  T-phase  amplitudes  can  be  modeled  approximately  using  a  simple  method  which  assumes  that  the 
amplitude  of  seismic/T-phase  coupling  at  any  location  is  dependent  upon  the  SOFAR  channel  depth, 
the  seawater  velocity  profile,  the  depth  of  the  seafloor,  and  the  distance  from  the  earthquake  epicenter. 
Future  plans  call  for  taking  into  account  propagation  effects  such  as  cylindrical  spreading,  blockage 
or  partial  blockage  of  acoustic  energy,  slope  effects,  and  seafloor  characteristics. 

Once  an  accurate  formulation  for  T-phase  amplitudes  has  been  determined,  it  will  be  possible  to  locate 

epicenters  by  correlating  observed  T-phase  coda  with  synthetics,  for  a  number  of  trial  source  locations. 

This  could  lead  to  an  accurate  method  of  source  location  using  only  two  receivers.  Note  that  it  is 

necessary  for  this  purpose  to  find  a  simple  expression  for  the  T-phase  coda  so  that  a  rapid  estimation  of 

source  location  may  be  made. 
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